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SUMMARY 


A decoupling control law synthesis technique is presented that integrates stability augmentation, decoupling, 
and the direct incorporation of desired handling qualities into an output feedback controller. The proposed design 
technique uses linear quadratic regulator (LQR) concepts in the framework of explicit model following (EMF) and 
proportional plus integral (PI) feedback. An idealized model is used as a part of the controller for robustness and 
handling qualities requirements. The PI error feedback is used to balance model-following performance with surface 
activity and to force steady-state errors to zero. The output feedback gains are then computed from a projection of 
the full-state gains. Closed-loop performance is shown by application of the control laws to the linearized equations 
of motion, and a six-degree-of-freedom simulation of an oblique-wing aircraft. Model following is shown to be quite 
good, but with some significant time delay in both linear and nonlinear evaluations. Decoupling of the longitudinal 
and lateral-directional axes is excellent, with some small oscillations evident in the nonlinear responses. 


NOMENCLATURE 


A,B,H,F 

On 

O y 

C 

DFBW 

EMF 

In 

i 

J 

K x 

K v 

LQR 

M 

OWRA 

PI 

P 

Q 

Q 

R 

r 

S 

u 

W 

x 


state space quadruple defining the model state and output equations 

normal acceleration, g 

lateral acceleration, g 

error selection matrix 

digital fly by wire 

explicit model following 

identity matrix of dimension n 

v/TT 

quadratic cost function 
state feedback gain matrix 
output feedback gain matrix 
linear quadratic regulator 
Mach number 

oblique wing research aircraft 
proportional plus integral 
roll rate, deg/sec 
quadratic state weighting matrix 
pitch rate, deg/sec 

quadratic control input weighting matrix 
yaw rate, deg/sec 

quadratic cross-weighting matrix for states and inputs 
control vector 

output vector weighting matrix 
state vector 


l. 


y 

output vector 

y e 

output error vector 

a 

angle of attack, deg 

0 

angle of sideslip, deg 

SaL 

left aileron deflection, deg 

&aR 

right aileron deflection, deg 

&hL 

left horizontal tail deflection, deg 

ShR 

right horizontal tail deflection, deg 

&R 

rudder deflection, deg 

$ 

roll angle, deg 


Subscripts 

e error 

Ie integral error 

m model 

p plant 

Superscripts 

T complex conjugate transpose 
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INTRODUCTION 

Early win d- t un nel and the oretica l studies have sho wn that a variable-skew oblique wing offers a substantial aero- 
dynamic advantage over conventional or symmetric variable-sweep aircraft for missions that requiroboth subsonic 
loiter and supersonic dash or cruise (Nelms, 1976 and Wilcr, 1985). The NASA Ames Research Center, Dryden 
Right Research Facility AD-1 flight program successfully demonstrated the oblique-wing concept and explored the 
low subsonic handling qualities and performance envelope (McMurtry, 1981). Because most of the aerodynamic 
advantage occurs at transonic and supersonic speeds, a high-speed flight research progra m was joi ntly sponsored by 
Ames-Drydcn and the U.S. Navy. The oblique wing research aircraft (OWRA) was to be based on a modified F-8 
digital fly-by- wire (DFBW) aircraft with a composite wing and wing-pivot mechanism substituted for the existing 
high wing structure (fig. 1). ■ . — - 

A major challenge in the OWRA program is the design and implementation of a control system architecture 
that will provide stabilization, decoupling, and acceptable fl ying qualitie s across the Mach number (A/), angle-of- 
attack, and wing-skew envelope. Because of the vehicle dynamic cross-coupling, the control system design must 
be applied to at least five (rigid body) degrees of freedom simultaneously, rather than separating the problem into 
the more typical two or three degrees of freedom in the longitudinal- and lateral-directional modes. The OWRA 
configuration provides an excellent opportunity to apply modern control methodologies to the problems associated 
with an oblique-wing configuration. 
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Model following has been shown to be an effective method for decoupling and stabilization of an initial OWRA 
configuration in previous work by Alag and others (1986). Implicit and explicit model-following (EMF) techniques 
were used successfully in a variety of controller structures to decouple the aircraft dynamics, but control surface 
activity was normally high, and steady-state model-following errors degraded the control system performance over 

a period of time. 

Vincent (1984) proposed an integrated model-following technique where dynamic elements are introduced into 
the constant gain controller to reduce steady-state model-following errors, and a nonoptimal output-feedback strategy 
is used to implement the resulting control law in a realizable form. In this way, handling quality requirements and 
trajectory tracking can be directly incorporated into a linear quadratic regulator (LQR) synthesis technique. Output 
error weighting is used instead of state weighting in the LQR cost functional, which gives the designer a better 
physical feel for the parameters being minimized. 

This paper presents a form of the integrated model-following technique as applied to the OWRA model. The 
resulting controller is structured as a proportional plus integral (PI) output-feedback control law. The effectiveness 
of the controller is shown by closed-loop time responses from linearized equations of motion as well as nonlinear 
six-degree-of-freedom simulation results at a given high subsonic flight condition and wing-skew position. 

PROBLEM DEFINITION 

Model-following is a useful technique when the dynamic characteristics for an ideal model can be specified. 
In EMF, the dynamic model is an integral part of the controller itself. The objective of EMF control is to force 
the aircraft to respond as the model would to a given pilot command. The model-following problem can be stated 

as follows: 

Given the linearized plant dynamic equation 

x p = ApXp + BpU p ( 

and the linearized model dynamics 

Xjji ■ 

where x p , x m , u p , and u„, are real vectors, and all matrices A p ,A m ,B p , and B m are of suitable dimensions, find the 
control Up that will force the aircraft states x p to approximate the model states x ra . Although sufficient conditions 
exist for guaranteeing “exact” or perfect model following, the conditions cannot be met by a physical system Ciyier, 
1964). This problem definition can be put in the context of the LQR, where the cost functional can be expressed in 
terms of the plant and model states or outputs. Therefore, the squared error between the plant and the model would 
be minimized, forcing the plant to follow the model states or outputs. An output error LQR formulation was used in 
this investigation primarily to simplify the task of selecting cost function weighting matrices and to give the designer 
better insight into the parameters being minimized. 

OUTPUT ERROR LINEAR QUADRATIC REGULATOR 

Consider the aircraft plant and model dynamic equations given above in equations (1), (2), and the out- 
put equations 


y p — HpXp + 


(3) 


y»n - + F m ll m 


(4) 


3 



The most common form of the quadratic cost is a function of the state and control vectors of the system, written 


as 



rOO 

J= x t Qx + u T i?u dt 

Jo 

(5) 

or in output formulation 

J = f y T Qy + u T i?u dt 

Jo 

(6) 


where Q > 0 and R > 0 are the quadratic weighting matrices. Define the output error as y e - y p — y m to determine 
the cost as a function of the error between the plant and model outputs. Integral error can be added by defining the 
error slate equation: 

*/ € = Cy c ( ? ) 

where C is a selection matrix operating on the output error vector. Thus the state and output equations can be 
augmented to form the following system: 


Xp 


1 

tu 
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O 
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Xp 


B p 

X/e 

= 
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X/e 
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( 8 ) 


y = 


y« 


yie 

= 

y m 



H v 

0 

0 


0 ~H n 

1 0 
0 H m 


H 



(9) 


Since the output equation (9) is a linear function of the states and controls, the cost functional equation (6) can then 
be restated as 

/ = 

with (^ > 0 and H > 0 , 5 > 0 , and given by 


r 

Jo 


x r Qx + x T Su 


u t 5 t x + u T flu dt 


( 10 ) 


Q = 


$ = 


R = 


H p J Q e H p 0 -HpQcHr, 

0 Qle 0 

-HrjQeHp 0 HjQ'Hm 

H p J Q e F p -Hp T QeF, 

0 0 
-H m T Q e F p HrrjQ'Frr 

Fp J Q e F p -F p T Q e F m ' 

L -F m J Q e Fp FrrjQ'Fn . 


+ R 


( 11 ) 


( 12 ) 


(13) 


where Q e is the weighting matrix for the output error terms, and Q u is the weighting matrix for the integral error 
terms. If (A, B) is stabilizable, and ( H, A) is observable, the linear steady-slate optimal control that minimizes J 
equation (10) is found by solving the steady-state Riccati equation (Bryson, 1975), and is of the form 


u = -K x \ 


(14) 


« 


* 
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The resulting optimal control law is a function of all the states in the augmented system. A projection from the 
state space to the output space is possible without changing the closed-loop eigenvalues (appendix) if the number of 
unique outputs is equal to or greater than the number of states in the augmented system. The resulting eigenvectors 
are not guaranteed to be unchanged, but the model-following control structure will dominate the response behavior 
by driving the error and integral error to zero. This projection can be stated as follows: 

Define a weighted output vector 

y = Wy = Hx + Fu 

where H = WH, and F = WF. 


Since H is assumed full rank, and (H T H) > 0, this implies that (H T H)~' exists. The output equation can 
then be rewritten in terms of the system states as 

x = (H T H)~ l H T y (15) 

The optimal control can therefore be written in the output form 

u = -[I-K x (H T H)- l H J Fr l K x (H T Hr l H i: y (16) 

V ■ * 

K, 


where K y can be partitioned as [ K p Ku K m ] (see fig. 2). 

It can be shown that if y = x (this implies H - W = I and F = 0), it is possible to reduce the above equation to 
the state form of the optimal control law equation (14). W and 12 can be chosen in such a way that only the desired 
output feedback paths will have nontrivial gain values. Thus, all the feedback loops to the model can be eliminated, 
preserving the desired closed-loop eigenvalues of the idealized model. The control system design procedure can be 
stated as follows: 

1. Determine idealized model dynamics, 

2. Choose quadratic cost matrices Q e ,R,Qie> 

3. Calculate optimal state feedback gains K x from steady-state Riccati equation, 

4. Choose output weighting matrix W so that only the desired feedback loops have nontrivial gains, 

5. Determine output feedback gains Ky from state gain projection, and 

6. Iterate steps two through five until desired system performance is achieved. 

Figure 2 shows a block diagram of the final PI feedback controller. The resulting control law is nonoptimal in a 
strict sense because of the gain projection, but the design technique incorporates some important features. 

1 . The output error formulation allows the designer to balance model-following performance with control surface 
activity, both in surface deflection and rate (if the actuator dynamics are modeled in the plant), 

2. The integral error feedback forces steady-state model-following errors to zero, thus making the control law 
more robust to modeling errors in the aircraft plant dynamics, 

3. The idealized dynamic feedforward model allows the incoiporation of lead-lag networks, stick prefilters, or 
general command shaping in the forward loop for handling qualities requirements, and 
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4. The augmented control system structure, a dynamic prcfilter with PI error feedback, is similar to classical PI 
controllers. Therefore, designers can apply knowledge and intuition from past design experience, an important 
factor in a highly interrelated problem. 

As in most methodologies, this synthesis technique has some drawbacks. 

1. The number of resulting gains, even in the output error formulation, is prohibitively high, and gain scheduling 
between flight conditions is a complex problem. 

2. Although the full state LQR has guaranteed phase and gain margins, the projection from state to output space 
negates any robustness guarantees except for closed-loop stability. 

3. The selection of the proper weighting matrices Q e and R is not a straightforward process. This problem is 
made more difficult because the desired time or frequency response behavior not specified in the model cannot 
be easily related to the weighting matrices. Additionally, great care must be taken when implementing a linear 
idealized model (perturbation model, valid only for small angle maneuvers, instantaneous response to inputs, 
etc.) in a nonlinear, six-dcgrce-of-frccdom flight system. 

It is hoped that future research can overcome these drawbacks and allow this technique to be applied to a broader 
class of problems. 

OBLIQUE WING RESEARCH AIRCRAFT DESCRIPTION 

The F-8 OWRA design configuration was based on a feasibility design study done by Rockwell International 
(White and others, 1984) which had its roots in the initial work done at NASA Ames Research Center in the early 
70s (Graham and others, 1985; Smith and others, 1975 and 1976). The F-8 OWRA is basically the F-8 DFBW 
aircraft with the 357 ft 2 wing replaced by a 200 ft 2 composite wing and pivot structure (fig. 1). Wing skew is right 
wing forward, from 0 (unskewed) to 65°. 

The OWRA aerodynamic model was based on a combination of wind-tunnel studies, computational fluid dy- 
namics, and the simulation aerodynamic database of the F-8 DFBW aircraft. Much of the data from all three sources 
had to be linearized, interpolated, or modified in some way to form a cohesive full-envelope aerodynamic model. 
This modified aerodynamic model was implemented in a six-dcgrcc-of-frecdom, nonlinear, fixed-base simulation 
used for piloted evaluation and control systems development. Linear models used in the control system design were 
obtained from this nonlinear aerodynamic model using the program LINEAR (Duke, 1987). A comparison of linear 
and nonlinear open-loop time histories showed excellent agreement indicating that the linear models could be used 
even at the maximum wing-skew conditions. 

The engine model for the F-8 OWRA is based on the simulation engine model of the F-8 DFBW. The engine 
is equipped with an afterburner and is modeled in the simulation as two main rotational elements (compressor and 
fan) co-Iocated along the thrust line in the vehicle x-z plane. 

Different types of trim exist because of the unique aerodynamics of oblique-wing configurations. Trim is de fined 
as the conditions necessary to set the derivative of the stale vector (x) to zero. When the wing is skewed, a nonzero 
sideslip angle 0 or roll angle 4> is required to keep the trim <j> = 0 = 0 simultaneously. For this reason, level trim 
is defined as wings level ($ = 0) varying 0 as trim r equir es, and s ideslip trim is defined as constant 0 (nominally 
zero) varying <£ as required. Figure 3 shows the relationship between trim 0 and trim 0 at a single flight condition. 

Although many potential trim strategics exist, the ailerons are used to trim the roll axis, symmetric horizontal 
tail is used to trim the pilch axis, rudder is used to trim the yaw axis, and the throttle is used to obtain a steady-state 
velocity. Inboard flaps are used for unskewed powered approach only. Each surface moves independently, allowing 
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the controls designer greater flexibility in surface management which becomes more important at the high skew 
conditions, where the ailerons are effective (but highly coupled) controllers in all three axes. 


OBLIQUE WING RESEARCH AIRCRAFT DESIGN EXAMPLE 


The OWRA design problem is well suited to the EMF design technique because a realistic model exists with the 
desired time response characteristics and closed-loop eigenvalues. The model used was a six-state approximation to 
the closed-loop dynamics of the OWRA at zero wing skew (table 2). A simple rate and acceleration feedback control 
system was designed for the unskewed aerodynamics and then implemented in a six-degree-of-freedom fixed-base 
simulation. The feedback gains for this simple control system were modified until satisfactory pilot comments were 
obtained. Linear models were then generated to approximate the closed-loop performance of this augmented system. 
The six states used were the minimum number required to incorporate four of the five classical rigid-body dynamic 
modes: short period, dutch roll, spiral, and roll mode. It was not necessary to include the open-loop phugoid mode 
in the desired model dynamics because it remains stable and relatively fixed throughout most of the flight envelope. 
The use of the zero-skew aircraft as a model has an additional benefit because the desired closed-loop dynamics 
of the aircraft remained constant over the wing-skew envelope. Thus, the pilot would not have to adapt to new 
characteristics as the wing is skewed. 

For the OWRA design example, the model and plant state vector is given by 


x = 


a 

0 

4> 

p 

Q 


r 


angle of attack , rad 
sideslip angle , rad 
roll angle, rad 
roll rate, rad /sec 
pitch rate, rad /sec 
yaw rate, rad /sec 


(17) 


and the control vector 


u 


ShL 

iR 
6 at 
&aR 


left horizontal tail deflection , rad 
right horizontal tail deflection, rad 
left aileron deflection , rad 
right aileron deflection , rad 
rudder deflection , rad 


(18) 


Actuators were not modeled to simplify the design process and to reduce the number of feedback gains required. 
Because of this, only surface deflection could be directly minimized in the design, as no cost was incurred because 
of the surface deflection rate. Numerical difficulties were encountered in the full-state feedback gain projection 
to output feedback when actuator dynamics were included in the augmented state vector. These difficulties were 
primarily a function of the limited number of output parameters that were selected for feedback, rather than a fun- 
damental limitation of the technique. The output parameters chosen were those typically used as feedback signals 
in classical control systems, especially those parameters related to the lateral-directional axis. Signal reliability and 
sensor redundancy were also factors in the selection. The output vector was identical for the model and the plant. 
The output error y e is given below as 


Pe 

9e 



One 
- °Ve 


roll rate error, rad /sec 
pitch rate error, rad /sec 
yaw rate error , rad / sec 
roll angle error, rad 
normal acceleration error, g 
lateral acceleration error, g 


(19) 
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with the error between the three angular rates of the plant and model 


yu = 


Pe 


r t 


roll rate error, rad /sec 
pitch rate error , rad /sec 
yaw rate error , rad /sec 


( 20 ) 


used as the input to the integral error feedback y Ie . The angular rates were chosen as integral error feedbacks to 
improve the steady-state model following and add additional damping to cross-axis perturbations. 


RESULTS 

The design flight condition was M = 0.8, altitude = 20,000 ft, and wing skew = 45°. This flight condition 
was chosen because it was contained in the proposed initial flight test envelope, and the cross-axis coupling was 
near a worst case for the entire flight envelope. Additionally, the total aircraft drag was near a minimum at the 45° 
wing-skew condition (fig. 4). Tables 1 and 2 show the system matrices for the plant and model respectively. 

Many design iterations were used to determine a set of weighting matrices that resulted in good system perfor- 
mance with low control surface activity. The final matrices used were the following: Q e = DIAG ([10,20,10, 500 , 
10, 100]), Qj e = DIAG ([500 ,500, 100]), and R = DIAG ([200 1 5 \ 10 7 / 5 ]). The large values in R are applied 
to u ro to eliminate any feedback loops to the model. The output weighting matrix W was constructed to eliminate 
all but the desired feedback gains. In this case, W = DIAG([ J 9 ; lO -5 !*])- The resulting output feedback gain 
matrices K p , Kj e , K m are shown in table 4. For this design case, it is interesting to note that the number of unique 
outputs is less than the number of states (Rank (H) = 9 , Rank (A) = 15) , yet the eigenvalues for the output 
feedback system remain close to the eigenvalues of the full-state system (table 3). 

The responses of the linearized open-loop aircraft model to a l°-ramp in elevator and aileron are shown in 
figures 5-6. The ramp begins at 1 sec and ends at 3 sec for each input. The open-loop aircraft responses are lightly 
damped (note a, q, and the accelerations) and the roll angle (<f>) due to pitch command is larger in magnitude than 
the roll angle due to roll command. The roll response in the pitch axis was found to be a strong function of angle 
of attack as well as nonlinear at this flight condition. For this reason, only small amplitude inputs were used in the 
design evaluations, since the linear models were not valid for large amplitude maneuvers. Time histories for the 
closed-loop system were run initially with linear aerodynamics and nonlinear actuator models with rate and surface 
limits included. 

After some design iterations, the control laws were implemented and tested in the nonlinear six-degree-of- 
freedom simulation. Figures 7-14 show the ideal linear model and closed-loop aircraft responses (both linear and 
nonlinear) to a 4°-roll ramp, while figures 15-22 show the responses to a 5°-pitch ramp. The closed-loop system is 
much more damped, and the roll due to pitch command has been all but eliminated. Good model following is evi- 
dent, but with some significant time delay evident for both linear and nonlinear cases; note <f> and lateral acceleration 
for the roll input, and q, a, for the pitch input. Oscillations can be seen in the nonlinear responses, especially in the 
cross-axis components, but the amplitudes are quite small. 


Ames Research Center 

Dry den Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, October 6, 1989 
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Table 1. Open-loop system matrices flight condition; M = 0 .8 , altitude = 20,000 ft, wing skew - 45°. 


A p 



' -0.7826 

0.0958 

0.0000 

0.0030 

0.9926 

-0.0003 


-0.0592 

-0.2908 

0.0387 

0.0259 

0.0001 

-0.9920 


0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0275 

- 

33.1431 

-53.6928 

0.0000 

-3.1250 

2 .0552 

1.7210 


-8.6816 

0.7975 

0.0000 

0.1679 

-1.0352 

0.1810 


-1.0092 

10 .7521 

0.0000 

-0.0213 

0.0080 

-0.7129 


Bp 


-0.0974 

-0.0974 

-0.0198 

-0.0302 

0.0000 

-0 .0166 

0.0166 

0.0008 

-0.0005 

0.0647 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

12 .9304 

-22 .2654 

15 .8467 

-11.8422 

13.2774 

-9.4073 

-10.8655 

-1.2311 

0.8797 

0.5694 

1 .9854 

-2.2579 

0 .5262 

-0.3276 

-6.2499 


Hp 


0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

20.1511 

-2 .4782 

0.0000 

0.0302 

0.1895 

0.0097 

-1.5285 

-7.4982 

0.0000 

-0.0426 

0.0036 

0.1972 


Fp 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

2.5097 

2.5097 

0.5103 

0.7802 

-0.0033 

-0.4291 

0.4291 

0.0216 

-0.0120 

1 .6674 
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Tabic 2. Desired model system matrices. 


-1.1568 

0.0000 

0.0000 

0.0000 

0.9566 

0.0000 

0.0000 

-0.2852 

0 .0387 

-0.0148 

0.0000 

-0.9600 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

-41 .7362 

0.0000 

-8.0000 

0.0000 

10.0000 

-16.8252 

0.0000 

0.0000 

0.0000 

-4 .0271 

0.0000 

0.0000 

10.9181 

0.0000 

-0.2726 

0.0000 

-3.3632 


-0.0848 

-0.0848 

-0.0494 

-0.0494 

0.0000 

-0.0166 

0.0166 

0.0000 

0.0000 

0.0647 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

11.7086 

-11.7086 

29.2528 

-29.2528 

8.1349 

-7.8228 

-7.8228 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-6.4014 


0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

28 .5286 

0.0000 

0.0000 

0.0000 

0.2342 

0.0000 

0.0000 

-7.3538 

0.0000 

-0.0394 

0.0000 

0.2075 


' 0.0000 0.0000 0.0000 0.0000 0.0000 ' 

0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 

2.1840 2.1840 1.2743 1.2743 0.0000 

-0.4291 0.4291 0.0000 0.0000 1.6674 . 
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Table 3. Closed -loop eigenvalues for full-state and output feedback. 


Full-state feedback 
-7.4002 ±4. 1541 * 
—6 .1927 ± 3.7359* 
-3.5624 ±2. 3180 » 
-0.2347 
-0.0341 


Output feedback 
-7.1095 ±4.8533* 
-5.3189 ± 4.1278* 
-3.3258 ±2.4188* 


-0.2271 

-0.0317 


- 0.0011 - 0.0011 


Model eigenvalues 

-2.5920 ±3.7464* 

Short period 

-3.5624 ±2.3180* 

Dutch-roll 

-0.0120 

Spiral 

-7.9137 

RoU 


Table 4. Model feedforward, error feedback, and integral error feedback gain matrices flight condition, 

M = 0 .8 , altitude = 20,000 ft, wing skew = 45°. 


K 


m 


K p 


0.4449 0.1051 

1.1031 -0.3790 

0.5542 0.2179 

0.5025 -0.1434 

0.7614 -0.5524 


0.1694 0.4585 

0.2076 0.3803 

0.1346 0.0891 

0 .0992 -0 .0900 

0.1042 0.0085 


0.0006 -0.0409 

0.0017 0.0613 

-0.0001 -0.0359 
-0.0021 0.0258 

0.0012 -0.0286 


-0.1858 -0.6725 

0.1436 0.7338 

-0.0243 -0.4903 
0.0101 0.3344 

0.5972 -0.1658 


0.0485 

0.0261 ' 

0.1025 

0.0051 

-0.0234 

0.0340 

0.0164 

-0.0275 

-0.1194 

0.1726 . 


0.0055 

-0.1032 ‘ 

0.0717 

0.0521 

-0.0200 

-0.0760 

-0.0179 

0.0476 

-0.0255 

-0.0145 


‘ -0.4697 

0.7455 

-0.1821 ' 

0.5597 

0.6492 

0.1385 

-0.3838 

0.1587 

-0.0805 

0.3108 

-0.3873 

0.0219 

-0 .4845 

0.2165 

0.3367 . 
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APPENDIX 


The following informal proof shows that the projection from state space to output space preserves the augmented 
system closed-loop eigenvalues if the number of unique outputs is greater than or equal to the number of states. 

Given the system equations 


x = Ax + B u (21) 

y = H\ + Fu (22) 

Where {A, B) is stabilizable, ( H, A) is observable, H is full rank, and Rank (H) > Rank (.4) . 

Given the optimal control in the form u = - K x x . Then the closed-loop system equation becomes the following 

x = (A-BK x )x (23) 

To determine the output form of the control, premultiply the output equation by ( H r H)~ X H J to get 


( H t H ) 1 H J y = x + (H r H) ' H T Fu (24) 

This inverse exists since H is full rank (by assumption). Substituting into the optimal control equation yields 

u = -K x (H r H)~ l H T y + K x ( H T H) 1 H T Fu (25) 

which can be rewritten in the final output form 

u = - [I - K X (H J H)- ] H J Fr l K X (H J H)~ l H 1 y (26) 

< v ' 

K, 

To examine the closed-loop eigenvalues of the output feedback system, substitute equation (26) into the system 
equations to yield 

x = Ax - BK y y (27) 

y = Hx - FK y y (28) 

Rewriting this last equation 

y = [I+ FK v ]~ l Hx (29) 

Thus the closed-loop output feedback system equation becomes 

x = (A- BK V [I+ FK v r l H)x (30) 

To show that the closed-loop eigenvalues of the state feedback (eq. 23) and output feedback (eq. 30) systems are 
the same, it is sufficient to show that the closed-loop dynamic equations of the two systems are identical. Comparing 
the two equations, it is necessary to show that K x = K V [I + FK y ]~ 1 H. 

From the definition of K v (eq. 26), 

[/ - K X (H T Hr 1 H r F)K v = Jr i (ff T tf)- , tf T (31) 
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( 32 ) 


Ky = K x { H r H ) -1 H t [ I + FKy] 
Collecting K y terms and post multiplying by H results in, 


K x = K y [ I + FKyV'H 


( 33 ) 
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Figure 3. Trim sideslip as a function of roll angle 
flight condition; M = 0.8, altitude = 20,000 ft, wing 
skew = 65°. 
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Figure 4, Total aircraft drag as a function of wing 
skew flight condition; M = 0.8, altitude = 20,000 ft, 
sideslip trim. 
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Figure 9, Closed-loop aircraft response to roll ramp 
sideslip angle. 



Figure 10. Closed-loop aircraft response to roll ramp 
yaw rate. 






Figure 13. Closed-loop aircraft response to roll ramp 
normal acceleration. 



Figure 14. Closed-loop aircraft response to roll ramp 
lateral acceleration. 
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Figure 15. Oosed-loop aircraft response to pitch 
ramp angle of attack. 
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Figure 16. Gosed-loop aircraft response to pitch 
ramp pitch rate. 
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Figure 17. Closed-loop aircraft response to pitch 
ramp sideslip angle. 



Figure 18. Closed-loop aircraft response to pitch 
ramp yaw rate. 




Figure 19. Gosed-loop aircraft response to pitch 
ramp roll angle. 



Figure 20. Closed-loop aircraft response to pitch 
ramp roll rate. 
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Figure 21. Gosed-loop aircraft response to pitch 
ramp normal acceleration. 
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Figure 22. Gosed-loop aircraft response to pitch 
ramp lateral acceleration. 
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